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We study long-lived resonances (lowest-loss modes) in hexagonally shaped dielectric resonators 
in order to gain insight into the physics of a class of microcrystal lasers. Numerical results on 
resonance positions and lifetimes, near-held intensity patterns, far-held emission patterns, and effects 
of rounding of corners are presented. Most features are explained by a semiclassical approximation 
based on pseudointegrable ray dynamics and boundary waves. The semiclassical model is also 
relevant for other microlasers of polygonal geometry. 
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I. INTRODUCTION 

A novel class of microlasers based on nanoporous 
molecular sieve host-guest systems has been fabricated 
recently by Vietze et al. . Organic dye guest molecules 
were put into the channel pores of a zeolitic microcrys- 
tal AIPO4 — 5 host. The aluminophosphate-crystals grow 
with natural hexagonal boundaries with a high degree of 
perfection. In terms of pump energy needed to reach las- 
ing threshold these microlasers can compete with semi- 
conductor based vertical-cavity surface-emitting lasers. 
This makes them a promising candidate for future appli- 
cations as e.g. optical communication devices. 

Microlasers and microresonators are not only relevant 
for experiments and applications but they are also of 
great interest from a theoretical point of view because 
they are open (or leaky) and they can be mesoscopic. In 
our case the latter depends on the lasing dye and on the 
crystal size. Both introduce a characteristic length scale: 
the wavelength A ranging from 600nm to 800nm, and the 
side length R of the hexagonal cross section of the crystal 
ranging from 2.6/im to 4.6/im in Refs. [|l], |^. In current 
experiments, larger crystals with R up to 20/im are under 
investigation. Only the ratio of these two length scales 
is relevant. It can be expressed by the dimensionless size 
parameter 20 < kR < 190, where k = 27r/A is the free- 
space wave number. For small kR the wavelength and the 
cavity size are of the same order. The system is micro- 
scopic. For large but finite kR the system is mesoscopic. 

The theoretical analysis of the lasing modes reduces 
to the two-dimensional problem of resonant modes in a 
passive dielectric microresonator of regular hexagonal ge- 
ometry (ignoring surface roughness). In Ref. 1^ prelimi- 
nary numerical computations of the full wave equations 
have been restricted to the near-field intensity pattern. 
We present a systematic numerical analysis of the hexag- 
onally shaped dielectric resonator with 20 < kR < 60, 
including resonance positions and lifetimes, near-field in- 
tensity patterns, far-field emission patterns, and effects 
of rounding of corners. To avoid convergence problems 



at corners, which had been a problem in Ref. M, we will 
employ the boundary element method (BEM) 

To study the deep mesoscopic regime kR > 60 we will 
introduce a semiclassical ray model. The semiclassical 
(short-wavelength) approximation is applied in the field 
of quantum chaos to relate quantum (wave) dynamics 
to their underlying classical (ray) dynamics. Most re- 
search efforts have been focused on closed resonators, so- 
called billiards, where the dielectric interface is replaced 
by a hard wall on which the wave function vanishes. The 
classical dynamics is free motion inside the billiard with 
specular reflections at the boundary. 

The regular hexagon belongs to the class of rational 
polygonal billiards. All angles (pj between sides are ra- 
tionally related to tt, i.e. = mjir/nj, where rrij, rij > 
are relatively prime integers. If all mj are equal to unity, 
for example in the case of the rectangle, then the dy- 
namics is integrable. In the presence of critical cor- 
ners with rrij > 1 the dynamics is not integrable but 
instead pseudointegrable Q. As for integrable systems 
the phase space is foliated by two-dimensional invariant 
surfaces [|[ ||. However, there are some peculiar fea- 
tures that distinguish these billiards strongly from inte- 
grable ones (and also from fully or partially chaotic bil- 
liards): (i) an invariant surface does not have the topol- 
ogy of a torus but instead that of a surface of higher 
genus H|, loosely speaking, a torus with additional han- 
dles; (ii) the dynamics is not quasiperiodic. This is re- 
flected, for example, by multifractal Fourier spectra of 
classical observables 0, Hi; (hi) the quantum-classical 
correspondence is exotic H; (iv) the quantum spectrum 
obeys critical statistics KQ]. To compute the quantum 
spectrum of a pseudointegrable billiard with a semiclas- 
sical treatment is extremely difficult, if possible at all. In 
the present paper, we will demonstrate that for a suffi- 
ciently open hexagonal dielectric resonator the spectrum 
and the eigenmodes can be computed semiclassically. As 
a byproduct, we gain an intuitive understanding of the 
numerical results. 

The paper is organized as follows. The system is de- 
fined in Sec. |l[ Section III presents the numerical analy- 
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sis. The semiclassical ray model is introduced in Sec. IV 
Finally, Sec. ^ contains our conclusions. 
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II. THE SYSTEM 

The numerical computation of lasing modes in a cav- 
ity with active medium and complex geometry is an ex- 
tremely difficult task; see Refs. |l^. Focusing on the 
effects of the geometry we here restrict our calculations 
to resonant modes in a passive cavity. Which of these 
resonant modes contribute to lasing is not important in 
our case since all long-lived resonant modes have very 
similar properties as we will see in the following. 

In the experiments on the microcrystals it has been 
shown that the electromagnetic field is TM polarized ||, 
^. Maxwell's equations reduce to a two-dimensional 
wave equation p3| 



(1) 



with wave number k and piece-wise constant index of re- 
fraction n(r). The index of refraction is n = 1.466 inside 
the cavity and 1 outside. The origin of the coordinate 
system r = {x, y) = (r cos 0, r sin 6) is located in the cen- 
ter of the hexagonal cavity. The complex-valued wave 
function tp represents the ^-component of the real-valued 
electric field vector Ez{r,t) — Ke[ip{r) exp {—iujt)] with 
= —1, angular frequency u — ck and speed of light 
in vacuum c. ip and V?/; are continuous at the boundary 
of the resonator. To model the situation in a laser, we 
impose the outgoing-wave condition 



^pr^h{0,k) 



exp (ikr) 



(2) 



for large r. With a real- valued n (passive resonator 
without absorption) this leads to modes that are expo- 
nentially decaying in time. The lifetime r of these so- 
called resonant modes or short resonances is given by 
r = — l/2clm(fc) with Im(fc) < 0. The lifetime r is re- 
lated to the quality factor Q — Re(ti;)r. The resonant 
modes are connected to the peak structure in scattering 
spectra; see e.g. Ref. Q. We are here only interested 
in long-lived resonances that provide a sufficiently long 
lifetime for the light to accumulate the gain required to 
overcome the lasing threshold. Note that extremely long- 
lived resonances are not relevant for lasers because they 
do not supply enough output power. 



A. Symmetry considerations 




FIG. 1: The eight symmetry classes of the hexagon. Odd (Even) 
symmetry is marked by solid (dashed) lines. 



way. Take two copies of a single H — a-mode. Rotate 
the first copy by 60° clockwise around the origin of the 
coordinate system and the second copy by the same angle 
counterclockwise. Subtracting the two gives a — -l-a-mode 
with the same k. The way from — h a to H — a, from 

+ + a to a, and from a to + + a is analogue. 

Hence, a-modes always appear in degenerate pairs having 
the same k. 

Each linear superposition of such a pair of degenerate 
a-modes, we denote them by ipi and ip2, is also a solution 
of the wave equation (|l]) . Because of the 6- fold symme- 
try of the system, we can find always two superpositions 
'ip± = "01 + pj-ip2 with the property that the correspond- 
ing intensities \ip+\'^ and £^re invariant under 60°- 
rotations. That means ip-^. and V'- change only by a phase 
factor expin± when rotated in real space by 60°. A full 
rotation about 360° does not change the wave function. 
Hence, n± = -Kq/S with q = -2, -1, 0, 1, 2, 3. Yet, q = 
and q — 3 are not allowed because the intensities of a- 
modes are not invariant under 60°-rotations. A rotation 
about 180° is identical to {x,y) (— x, — y). From this 
we find Cl± — ±7r/3 for 'ip± formed by a-modes of type 
H — , — h and Q± — ±27r/3 for ij;± formed by a-modes of 
type H — h, . These requirements lead to 



P± 



exp (z»±)0i(r) - V;i(f) 
exp (if]±)V'2(r) - V'2(f) 



(3) 



where r ^ (0, 0) is an arbitrary point and r is r rotated 
by 60° counterclockwise. Formula (||) is valid for any 
relative phase of ipi and 'ip2 ■ 

The 6-modes cannot be converted into each other since 
each 6-mode is invariant under 60°-rotations up to a 
phase. Hence, 6-modes do not form degenerate pairs. 



Figure |l| shows the eight symmetry classes of the 
hexagon; see e.g. Ref. |l^. In the notation — ha, . . ., the 
first sign is + if the wave function is even with respect to 
X —^ —X, and — otherwise. Correspondingly, the second 
sign refers to y —^ —y. The letter a indicates two sym- 
metry lines, whereas the letter b indicates six symmetry 
lines. 

The H — a and the — h a-modes have exactly the same 
complex k. The reason is that a H — a-mode can be 
converted to a — (- a-mode and vice versa in the following 



B. Rounding 

For our numerics it will be necessary to round the cor- 
ners of the resonator slightly as depicted in Fig. ^. In 
terms of polar coordinates the parametrization of the 
boundary reads 



(cos 6* - ^ sinOy + (^ sinOy + {cos 9+ smOy 

(4) 
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where s is the rounding parameter, a positive even in- 
teger, s = 2 gives the circle, whereas s — > oo gives 
the hexagon with flat sides and sharp corners. The 
parametrization in Eq. preserves the full symmetry 
of the problem. Note that the parametrization used in 
Rcf. § is different. 




FIG. 2: Shape and discretization (circles) of the resonator with 
(a) s = 20, (b) s = 200 and 7] = 0.1. For iUustrational reasons only 
A'^ = 200 discretization points are used. 

It is easy to show that the maximum curvature for 
the rounded hexagon in Eq. (Q) is given by (s — l)/3i? 
for s > 6. The radius of curvature p per wavelength 
A = 27r/Re(A:) is therefore given by 

p 3 Re{kR) 

for large s. 

III. NUMERICAL ANALYSIS 

A frequently used numerical method to solve wave 
equation (|l|) with outgoing-wave condition (|^) is the 
wave- or mode-matching method 0] . The wave function 
is expanded in integer Bessel functions inside the cav- 
ity and in Hankel functions of first kind outside, so that 
the outgoing-wave condition (§) is fulfilled automatically. 
The Rayleigh hypothesis asserts that such an expansion 
is always possible. However, it can fail for geometries 
which are not sufficiently weak deformations of a circu- 
lar cavity [l^ . So, using the mode-matching method 
for the hexagonal resonator in Ref. |^ is questionable. 
And indeed, according to Rcf. serious convergence 
problems appear in this situation. To avoid this problem 
we apply the BEM Exploiting Green's identity, the 
two-dimensional differential equation is replaced by a 
one-dimensional integral equation defined only along the 
boundary of the cavity. To be slightly more precise, the 
boundary appears twice in the formalism: one bounds 
the interior of the cavity and one bounds the exterior. 
The boundaries are then discretized by dividing them 
into 2N small boundary elements. For the boundary (^) 
we choose the following discretization 

9 = ^ + rism[6{^ + 7T/2)] . (6) 

The discretization points on the circle ^ G [0, 2tt) are cho- 
sen to be spaced equidistantly. The parameter rj deter- 
mines the distribution of points on the cavity (the shape 



of which is completely determined by the rounding pa- 
rameter s). 7] = gives an uniform density of points 
while 77 > gives an enhancement near the corners as 
illustrated in Fig. |[ We always use 77 = 0.1. 

A. Resonance positions in the complex plane 

Figure || shows the total cross section a as function 
of the dimensionless wave number kR [Im(fci?) = 0] 
for plane- wave scattering with incidence angle 9 — 15°. 
For numerical details of how to compute this quantity 
see Ref. We observe a single-mode spectrum, a 

series of roughly equidistant peaks (the fine structure 
around kR « 24.5 is an artifact of the BEM which can 
be removed by increasing the number of discretisation 
points ) . We will count the peaks by the mode index (or 
quantum number) m according to increasing kR starting 
with m = at fci? = 0. The first peak at kR « 20.5 
in Fig. H is then labeled by to = 23. The position of a 
given peak kR and its width F in Fig. |3| are related to the 
complex value of kR^odc of the corresponding resonant 
mode by means of fci?,„odc ~ kR — ir/2. The subscript 
"mode" will be dropped from now on. 
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FIG. 3: Calculated total scattering cross section cr/R vs. kR for 
a hexagonal resonator. The plane wave is incidence at 15° to the 
horizontal side faces, s = 100, and 2N = 2000. 

The BEM not only solves the scattering problem but 
it also can compute the resonant modes. In Fig. ^ we 
plot the real and imaginary part of kR of a resonance as 
function of the rounding parameter s. Re(fci?) saturates 
around s = 55, whereas Im(fci?) saturates at s > 100. We 
take s = 100 translating to p/A w 0.11. It is surprising 
that one has to decrease p to values one order of mag- 
nitude smaller than the wavelength. In the following we 
will fix p/A w 0.11, that means when we change Re(fci?), 
we have to change s accordingly. 

Figure ^ shows the long-lived resonances in the com- 
plex plane inside the strip 20 < Re(fci?) < 60. Our 
numerics cannot cover the full regime of the microlaser 
experiments |], | 20 < Re(fci?) < 190. Small values 
of |Im(fci?)| correspond to long-lived modes and, corre- 
spondingly, large values correspond to short-lived modes. 
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FIG. 4: Real and imaginary part of kR as function of the rounding 
parameter s. 2N = 2000, m = 26, kR = 22.8725 - j0.1064 for 
s = 100. 

Four features in Fig. || are striking: (i) the a-modes 
are two-fold degenerated as predicted from the sym- 
metry considerations, (ii) The 6- modes come in quasi- 
degenerated pairs with shghtly different kR (differences 
in Re(fci?) and lm(fci?) are of the same order which can- 
not be seen in Fig. H). It is not possible to resolve 
these small splittings in the scattering cross section in 
Fig. |. (iii) The values of Re(fci?) are approximately 
equidistantly spaced, (iv) Highly-excited resonances with 
Re(fci?) > 34 lie approximately on a smooth curve. 
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FIG. 6: (a) Wave number Re(feiJ) vs. mode index m. The data 
can be approximated by a straight Hne Re(fci?) = vi^m -\- mo) with 
V Si 0.8238 and mo ^ 1.7516. (b) Im{fci?) Re{A:i?) vs. m. For 
m > 34, the data is well approximated by the mean value —2.87. 

From Fig. ^(b) we can infer that the product of real 
and imaginary part of fci? for highly-excited modes does 
not depend on the mode index m. We find 



Im(fci?)Re(fci?) « -2.87 



(7) 



for n = 1.466. Later we will see that this remarkable re- 
lation is a good approximation not only in the numerical 
accessible regime but also in the full experimental regime. 
Relation (|^) implies that the lifetime is proportional to 
Re(fc)i?^ and the quality Q is proportional to Re(fci?)'^ 
(Q ranges from roughly 70 to 6400). This finding is rel- 
evant for the future experiments on microcrystal lasers: 
it implies that the laser threshold decreases as the size of 
the resonator is increased. 

Due to the qualitative change in behaviour around 
Re(fci?) K. 34 in Figs. | and |(b) we distinguish the re- 
gions Re(fci?) < 34 and Re(fc_R) > 34. We refer to the 
former one as the microscopic regime and to the latter as 
the mesoscopic regime. 



FIG. 5; Long-lived resonances in the complex plane, cf. Fig. |l[ 
Filled circles on the solid line are the semiclassical solutions in 
Eqs. (|l^) and Q from Sec. [v[ 

The equidistant spacing of Re(fci?) can be seen more 
clearly in Fig. ^(a). Pairs of (quasi-) degenerated modes 
are labeled by the mode index. We see that all pairs 
lie extremely close to a line Re(A:i?) = v(rn -I- mp) with 
mode spacing v and shift mo- Linear regression gives v « 
0.8238 and mp ~ 1.7516. The same mode spacing can be 
observed in Fig. 0. The corresponding free spectral range 



B. Mode structure 

Figures |^(a) and (b) show the near-field intensity pat- 
tern of the resonances 26 -I- +a and 26 a (mode index 

m = 26), respectively. While these two resonances are 
standing waves (ignoring the uniform temporal decay), 
the corresponding superpositions ^+ and ^_ in Figs, ^(c) 
and (d) are unidirectional traveling waves. We call them 
chiral modes. Such a mode is specified by its mode in- 
dex m and a label -I- (traveling counterclockwise) or — 
(traveling clockwise). In the following we will deal only 
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with these chiral modes which are more straightforward 
to compare to the ray dynamics. We remark that the 
Husimi function representation frequently appUed to op- 
tical microcavities (see e.g. Ref. [^) does not provide 
more insight in our case. 




FIG. 7: Calculated near-field intensity pattern |?/)(r)p of two-fold 
degenerated resonances, (a) 26 + +a and (b) 26 a with well- 
defined parity, (c) 26-1- and (d) 26— with well-defined chirality. 
Intensity is higher for redder colours, and vanishes in the dark 
regions. kR = 22.8725 - i0.1064, s = 100, 2N = 2000. 

Figure || shows a higher-excited superposition . The 
following properties can be observed: (i) the intensity is 
concentrated along the boundary of the cavity, resem- 
bling whispering- gallery modes in circular or weakly de- 
formed circular cavities, (ii) The wave pattern looks reg- 
ular. An approximate nodal-line structure with a pecu- 
liar twist is visible, (iii) The emission is predominantly 
at the corners. This is in agreement with the laser emis- 
sion measured in the experiments in Refs. ^ |2|. (iv) 
Outside the cavity the light propagates along certain di- 
rections not parallel to the facets. The latter fact can 
be seen more clearly in the far-field distribution shown 
in Fig. There are six emission peaks with angular 
width of « 14° and angular distance to the nearest facet 
of « 17°. 

The sensitivity to rounding found for the resonance po- 
sitions in the complex plane, see Fig. also shows up in 
the mode structure. Figure |9| shows the near-field inten- 
sity pattern of a rounded hexagon with s — 20, i.e. p « A. 
The emission is again at the corners, but it is reduced. 
Moreover, the directionality has decreased and the high- 



FIG. 8: Chiral resonance 50- 
2N = 4000. 



kR = 42.6318 - j0.06766, s = 200, 



intensity directions are now parallel to the edges. The 
latter fact can be seen better in the far-field emission 



pattern in Fig. 11 



o 



FIG 
Fig.E 



9: Resonance 50— in a rounded hexagon with 
kR = 42.7099 - j0.01836, 2N = 4000. 



20, cf. 
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FIG. 10: Far-field emission pattern, for large r, of (from 

above) 50 a, 50 + +a, 50—, and 50+ modes; cf. Fig. H. Vertical 

lines mark the directions parallel to the edges. 
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A. Geometric optics 

Geometric optics follows from wave equation in the 
limit Re(fci?) oo neglecting all interference effects. 
In the following we focus on low-index materials with 
"-mill < n < n,„ax where = 1/ sin 60° w 1.16 and 

"max — 1/ sin 30° = 2. The lower bound guarantees 
that a six-bounce periodic ray with angle of incidence 
di = 60° > 9c (with sin 6'c = 1/") is trapped within the 
hexagonal resonator by total internal reflection at the 
facets; see Figs. |l^ (some aspects of the much simpler 
case n < Umin have been studied in Ref. Q). The pe- 
riodic ray is marginally stable with respect to shifting 
it along the boundary. In this way we obtain a whole 
family of periodic rays with identical length and angle of 
incidence. 

The upper bound n^ax — 2 ensures that triangular- 
shaped periodic rays with 30° angle of incidence are not 
totally reflected. The periodic-ray family with 9i = 60° is 
then the only long-lived family of short period. Yet, since 
periodic rays are dense in phase space, we cannot exclude 
the possibility that there are also long-lived periodic rays 
of high period. 

It should be emphasized at this point that the (infi- 
nite) long lifetime of the rays is relevant for our purpose, 
whereas the periodicity of the rays is not relevant. The 
latter comes in here simply because periodic rays (and 
their neighborhood) can be long-lived. We will not apply 
methods from semiclassical periodic-orbit theory |23[ | . 

The periodic rays are unstable with respect of chang- 
ing the angle of incidence. 



Figure 
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shows a ray with 
slightly different initial angle of incidence. The ray is 
slowly diverging from the central one. After some time, 
it reaches the corner on its other side at (almost) nor- 
mal incidence. Consequently, it then escapes refractively 
(with probability close to 1). 




FIG. 11: Far-field emission pattern of the resonance (from above] 
50 a, 50-f +a, 50—, and 50+ in the rounded hexagon; cf. Fig. bl 



IV. SEMICLASSICAL ANALYSIS 



Having computed numerically the long-lived reso- 
nances in the regime 20 < Re(fci?) < 60, we now in- 
troduce a semiclassical ray model in order to see what 
happens in the deep mesoscopic regime Re(fci?) > 60. So 
far, semiclassical approximations of dielectric resonators 
have been discussed only for the case of smooth bound- 
aries ||l^, |2^, |l] . Our heuristic approach is divided into 
three steps: geometric optics, semiclassical quantization 
and emission mechanisms. 




FIG. 12: Semiclassical ray model. Thick line marks a member of 
the family of long-lived rays, other members are obtained by shift- 
ing the ray along the boundary (not shown). The thin line marks 
a ray with slightly different angle of incidence. Arrows indicate 
emission due to pseudointegrable dynamics (thin) and boundary 
waves (dashed). 

As long as the ray does not escape from the cavity 
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it behaves as if it were in a closed resonator, i.e. in a 
billiard. It is illuminating to examine the ray dynamics 
of the open system in terms of the invariant surface of the 
hexagonal billiard. The genus of the invariant surfaces is 
given by the general formula [Q 



9 2^ ' 



- 1 



5=1 



where the sum is over all corners with angles 



(8) 



jT^/rij and Af is the least common multiple of the n^. Ain 



For the hexagon rrij = 2 and rij = 3 for all j; hence, 
J\f = 3 and finally g = 4. The surface of genus four 
is shown schematically in Fig. |l3|. The topology is the 
same for all initial conditions. First, we consider the spe- 
cial case where the surface is foliated by periodic orbits. 
Each periodic ray with 9i = 60° appears here as two 
disjoint circles with different sense of rotation. These 
two circles correspond to the two unidirectional travel- 
ing waves. In the following we will focus on the circle 
with clockwise rotational sense. On the same invariant 
surface we find also periodic rays of "bouncing-ball" type 
with angle of incidence 9i = 0°. In the open system, the 
bouncing-ball rays are not long-lived since they are not 
totally reflected. 

The dynamics on a generic invariant surface in the 
neighborhood of the special surface discussed above is 
as follows. A nonperiodic ray starting near the central 
periodic ray with slightly different initial angle of inci- 
dence stays in the vicinity of the central ray by winding 
around the handle many times. Finally, it has to leave 
the vicinity of the periodic ray since pseudointegrable 
motion is ergodic on generic invariant surfaces p4[ . Be- 
fore the ray can reach another handle (in the grey region 
in Fig. and perform complicated dynamics, it escapes 
refractively from the cavity. Hence, the sufficiently open 
resonator does not see the complicated topology of the 
full invariant surface but just two disjoint tori. In this 
sense, the openness moves the system closer to integrabil- 
ity. That is the reason why we can derive in the following 
subsection a practical semiclassical approximation for the 
open hexagonal resonator whereas this is impossible for 
the hexagonal billiard. 




FIG. 13: Invariant surface of the hexagonal billiard. The small 
circle on the_lcft (right) hand side represents the central periodic 
ray in Fig. cycling (counter-) clockwise. In the grey region the 
condition for total internal reflection is not fulfilled. 



B. Semiclassical quantization 

Having identified the long-lived part of classical phase 
space, we now quantize it in a semiclassical approxima- 
tion. This approximation is valid for small 1/Re(fci?) 
and it does include interference effects, in contrast to the 
geometric-optics limit. 

Because of the absence of stability of the periodic rays 
the conventional approach in laser optics, the paraxial 
approximation (see e.g. Ref. (2^), does not work. Never- 
theless, the semiclassical quantization is simple. The idea 
is that an integer number m = 1, 2, 3, ... of wavelengths 



X/n = 2n/nRe{k) fits on the path length of the 



long-lived ray shown in Fig. |1^, taking into account the 
phase shifts at the dielectric boundary. It will turn out 
that the number m defined in this way is identical to the 
mode index m as used in Sec. [II. 



To compute the phase shifts we consider the text-book 
treatment of reflection of a plane wave at an inflnitcly 
extended dielectric interface; see e.g. Ref. This 
simplified setup is justified if A ^ i?, i.e. Re(fci?) ^ 1. 
We shift the origin of the coordinate system such that 
?/ = is the dielectric boundary bounding the lower- 
index region from below and the higher-index region from 
above. In the higher-index region there is an incident 
wave 

ipi = Aiexp[iRe{k)n{x sin 9i + y cos 9i)] , (9) 

with amplitude Ai. We set Ai — 1 without loss of gener- 
ality. The reflected wave is given by 

ipr = Ar exp [iKe{k)n{x sin 9r — y cos 9r)] ■ (10) 

In the lower-index region there is an emitted wave 

= Aeexp[iRe{k){xsm9e + ycos9e)] ■ (11) 

The boundary conditions 

Mx,0)+Mx,0)^Mx,0) (12) 

and 



dy 









dy 

{x,0) 




dy 



lead to 9i — 6*.^, Snell's law n sin 9i — sin 9^, Ai + A,. 
and the Fresnel formula 



A-rp 



1 — ia 
1 + ia 



with 



\/ sin^ ' 



- 1 



(13) 



A, 



(14) 



(15) 



n cos 9i 

The quantization (or resonance) condition then reads 
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with the length of the periodic rays I = 3^/3R. 
some algebraic manipulations we arrive at 



Re{kR) 



2tt 



with the total boundary phase shift /? given by 



tan —P = a 
6 



After 



(17) 



(18) 



(17|) explains the single- 
5(a). With n = 1.466 the 



The quantization condition 
mode spectrum in Figs, and 

mode spacing is 27r/3A/3n « 0.8248 in agreement with 
1/ = 0.8238 obtained by fitting to the numerical data. 
However, the shift (3 « 1.562 differs a bit from the nu- 
merically obtained shift toq ~ 1.7516. To understand 
this discrepancy we plot in Fig. th e phase shifts f3„i 
defined as (3 computed from Eg. (|17|) inserting the nu- 
merically computed values of Re(fc^. Clearly, going to- 
wards the semiclassical limit m ^ oo leaves the crude 
fitting value ttiq « 1.7516 in favour of a smaller value 
/3oo = linim^oo Pm closer to our semiclassical prediction 
/3« 1.562. 
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FIG. 14: m VS. mode index m. 
sical prediction fj 1.562. 



Dashed line marks the semiclas- 



Let US assume that /?,„ has the following form 



7 m 



(19) 



and that /3oo equals our semiclassical solution /3 k, 1.562. 
Figure |l^ shows A/3 = f3m — Poo as function of m > 34 in 
a log-log plot for fixed n = 1.466. Linear regression gives 
7 w 0.89 and S « 0.5038. Based on this numerical finding 
we conjecture that the next order in the semiclassical 
approximation (|l^) is of the form j/y/m. 

Observation of Figs. |^(c) and (d) shows that the "nodal 
lines" close to the boundary of the cavity are not ex- 
actly parallel to the boundary but slightly tilted (about 
5°). This angular shift decreases as one goes to higher- 
excited modes. For example, it is not visible by eye in 
Fig. 1^. Specular reflection of a plane wave at the bound- 
ary leads always to nodal lines parallel to the boundary. 



-1.8 p 
-1.9 
-2 

S 

< -2.1 
-2.2 
-2.3 
-2.4 




1.2 1.25 1.3 1.35 1.4 1.45 1.5 
ln(m) 

FIG. 15: In A/3 vs. Inm writh m > 34. For quasi-degenerate modes 
the arithmetic mean of A/3 is taken. 



Hence, it is reasonable to interpret the angular shift as 
a kind of non-specular reflection. The angular shift can 
be traced back to the fact that the wave function re- 
stricted to the boundary is a periodic function. Because 
of the periodicity, the Fourier spectrum of the wave func- 
tion along the boundary is discrete. The discrete peaks 
can be related to a discrete set of allowed angles of in- 
cidence. The angle 9i — 60° typically falls between two 
neighboring discrete angles. In Sec. |A| in the Appendix 
we compute the angular shift analytically to first order 
as AO = 9i — 60° = 0^ — 7r/3 « — \/3/3/r7i in agreement 
with the numerics (not shown). In the semiclassical limit 
TO — > 00 the angular shift vanishes. Our angular shift 
is different from those cases discussed in the context of 
non-specular refiections (see e.g. Ref. |^), in that it 
happens above the critical angle of total internal reflec- 
tion. Moreover, it has nothing to do with the angular 
Goos-Hanchen effect (see e.g. Ref. |2^]), which is the 
Goos-Hanchen shift at a curved interface. 

The number of wavelengths fitting on a path length 
determines the symmetry class. Straightforward consid- 
erations give Tab. |[ The semiclassical quantization pro- 
cedure neither distinguishes between -| — and — -|- par- 
ties nor between -|--|- and parties. Hence, within the 

semiclassical approximation these modes are exactly de- 
generated, regardless whether they are of type a or 6. 
This explains the numerical finding of quasi-degenerate 
6-modes. 

The fact that our model predicts a single-mode spec- 
trum, a one-parameter to = 1,2,... family of long-lived 
resonances, is related to the fact that there is only one 
long-lived family of periodic rays (with short period). 
There is no "transversal mode index" which would be 
expected if the paraxial approximation were applicable. 
Higher transverse modes correspond to rays with angle 
of incidence considerably different from 60°. They are 
much shorter-lived which becomes clear in the following 
subsection when we discuss the emission mechanisms. 
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m 


symmetry class 


1 


H a, ha 


2 


+ + a, a 


3 


+ -b,- + b 


4 


+ + a, a 


5 


H — a, — ha 


6 


+ + b,--b 



we estimate the outgoing relative intensity due to pseu- 
dointegrable leakage for the hexagon as 



AT 47rc P{n) 

P 3Re(fc)i?2 n2 



(22) 



TABLE I: Symmetry class and mode index m. The table is peri- 
odic in m with period 6. 



C. Emission mechanisms 

Having derived a semiclassical quantization condition 
for the real part of kR, we now compute its imaginary 
part. The temporal behavior of the intensity of a reso- 
nant mode is / oc exp [2 Im(ct;)t]. The outgoing relative 
intensity per unit time is 



The derivation can be easily extended to any regular 
polygon. Cavities with integrable internal dynamics, like 
the equilateral triangle and the square, have A/p ~ 0. 
That means the neighborhood of the long-lived family of 
periodic rays has roughly the same lifetime as the peri- 
odic rays. In such a case, we expect a multimode spec- 
trum in the mesoscopic regime. Indeed, this has been 
found in scattering experiments on the dielectric square 
in Ref. |2|. 

Remarkably, both contributions in Eqs. (|2|) and ( [2^ ) 
have the same dependence on Re(fc). In either case, the 
contribution vanishes in the limit Re(/c) — > oo, reflecting 
their wave nature. Adding both contributions. 



AI = A/i 



bw 



AI 



1 dl 
l~di 



= — 2Im((jj) = — 2clm(fc) . (20) gives the central result 



What are the mechanisms for this decay of intensity? 
Obviously, there is no classical mechanism. In the frame- 
work of ray optics, the periodic rays in Fig. |l^ remain 
forever in the cavity. Hence, we have to include explic- 
itly wave effects. We identify three candidates of such 
wave effects, all of which are related to the corners. The 
first one is diffraction at corners. Corner diffraction may 
be a emission mechanism since a wave with finite wave- 
length coming close to a corner is diffracted partly to the 
exterior and partly back into the interior with another 
spectrum of directions for which the condition of total 
internal refiection may not be fulfilled at the next reflec- 
tion. We will argue that corner diffraction is important 
for the emission directionality but not for the escape rate. 
For the escape rate we find two new relevant effects which 
are responsible for transport to corners. We call them 
boundary-wave leakage and pseudointegrable leakage. 

The boundary-wave leakage is illustrated in Fig. 
An evanescent boundary wave travels along an infinitely 
extended dielectric interface from — oo to -hoo. We as- 
sume that at a finite interface the boundary wave fully 
separates from the interface at the corner. In Sec. ^ in 
the Appendix we determine the outgoing relative inten- 
sity for any regular polygon (equilateral triangle, square, 
hexagon, ...) due to boundary waves. In particular, for 
the hexagon we find 



Im(fci?) Re(fci?) = /(n) 



with 



f{n) = fhw{n) + fp{n) , 



fhwin) 



3n3 



8\/3n74- l(n2 - 1) 



and 



27r/3(n) 



(23) 



(24) 



(25) 



(26) 



(27) 



A/bw = 



3c 



4Re(fc)i?2 ^3^2/4 _ 1(^2 _ 1) 



(21) 



The pseudointegrable leakage is due to the fact that 
wave optics does not realize exactly 9i = 60° as already 
mentioned and shown in Sec. ^ in the Appendix. Putting 
this small angular deviation A0 into the initial conditions 
gives rise to rays with finite lifetime due to the pseudoin- 
tegrable dynamics; see Fig. |l^. In Sec. |^ in the Appendix 



For n = 1.466 we find lm{kR)Re{kR) w -2.837 in rea- 
sonable agreement with the numerical data in Figs. |^ 
and ||(b). Figure [l^ compares the semiclassical result 
in Eqs. (|2^)-([27|) to the numerical data as function of 
n. While some resonances are very well described by 
the semiclassical approximation, some are only roughly 
described. The latter cases correspond to the strong fluc- 
tuations around the semiclassical hyperbola (|2j) already 
visible in Fig. ^j. 

Our ray model also explains some aspects of the mode 
structure in Figs. || and |l^. The whispering-gallery-like 
structure is related to the geometry of long-lived rays; see 
Fig. ^ The envelopes of the emission direction are de- 
scribed by the boundary-wave leakage and the pseudoin- 
tegrable leakage; cf. Figs. ||, ^ and ]l2|. The distribution 
within the envelopes is not predicted by our model. Here, 
corner diffraction is important because both the bound- 
ary waves and the pseudointegrable ray dynamics lead to 
escape arbitrarily close to the corner. 

Our ray model gives a natural explanation for the 
sensitivity to rounding of corners. Let us assume that 
a sufficiently small rounding only slightly weakens the 
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DC 

DC 
DC 




an intuitive way, except the emission directionality. To 
describe the directionahty properly it is necessary to con- 
sider corner diffraction in the future. 

Our results are not only relevant for microcrystal lasers 
but also for other kinds of polygonal-shaped microlasers. 
To name a few: hexagonal lasers with attached op- 
tical waveguides based on semiconductor heterostruc- 
tures [^0|, hexagonally shaped solid polymer dye mi- 
equilateral-triangular laser cavities fab- 
semiconductor heterostructures [B2|, and 



31 



crocavities 
ricated from 

square laser micropillar cavities based on dye-doped poly- 
mers [BSl. 



FIG. 16: Im(kR) Re(fciJ) vs. n for several resonances with m > 50 
(solid curves). The contribution of the boundary waves, Eq. (|2£|), is 
shown as dotted curve The entire contribution, i.e. with pseudoin- 
tegrable leakage, Eq. (ETp, is shown as dashed curve. The vertical 
line marks n = 1.466. 



loss due to boundary waves, whereas it reduces corner 
diffraction strongly. Numerical simulations of ray dy- 
namics in rounded hexagons show that the periodic rays, 
see Fig. ^ are stabilized. This statement is indepen- 
dent on the particular chosen boundary parametrization. 
Hence, there is no pseudointegrable leakage. This leads 
to considerably larger lifetimes and emission parallel to 
the facets; compare Figs. ^^ll[ 



V. CONCLUSION 

We have discussed the properties of long-lived reso- 
nances in hexagonally shaped dielectric microcavities. 
These microcavities play an important role for experi- 
ments on a class of microcrystal lasers ||l|, ^ . 

The numerical analysis revealed the following facts: (i) 
The resonance positions in the complex plane are approx- 
imated by a hyperbola. This fact is relevant for future ex- 
periments on microcrystal lasers because it implies that 
the laser threshold is lower for larger crystals, (ii) The 
near-field intensity pattern show a whispering-gallery-like 
structure, (iii) The emission is at corners in agreement 
with the experiments, (iv) The emission is directed. The 
high-intensity directions are not parallel to the facets, (v) 
The lifetimes and the emission directionality are sensitive 
to rounding of the corners. 

The numerical analysis does not cover the full exper- 
imental regime. To overcome this limitation, we have 
introduced a semiclassical approximation which can be 
easily extended to any cavity of regular polygonal ge- 
ometry. Our semiclassical ray model contains two new 
emission mechanisms: leakage due to boundary waves 
and due to the pseudointegrable ray dynamics. Explicit 
expressions for the resonance positions can be given even 
though the system is not integrable. The semiclassical 
approximation explains all numerical findings (i)-(v) in 
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APPENDIX A: ANGULAR SHIFT 

In this section we discuss the semiclassical deviation of 
the angle of incidence from 60°. Let us first investigate 

and H — h- modes. Along the boundary, these modes 

are periodic with period 3i?, i.e. •ip{s + 3i?) = tpi-s), 
where the arclength s € [0, 6i?) parametrizes the bound- 
ary. Hence, the wave function along the boundary can 
be expanded as the following Fourier series 



E 

Af=-oo 



Am exp 



, , 
I — Ms 
3i? 



(Al) 



From this we see that the conjugate momentum to s or 
in other words the momentum component parallel to the 
boundary is "quantized" according to 



iRe(fc) ! 



= — M 
3i? 



(A2) 



with integer M. Inserting the quantization condition (|ri 
leads to 



sm ( 



= V3- 



M 



Linearising this equation around 9 ^ 9i 
yield in the semiclassical regime m ^ /3 



60° 



m m 



(A3) 
= 7r/3 

(A4) 



The smallest Ad is realized for M = m/2 (m is an even 
integer for and H— I— modes, see Tab. |), 



AO 



3nRe{kR) 

Analogous arguments concerning — h and 
give the same result as in Eq. (A5). 



(A5) 
-modes 
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APPENDIX B: BOUNDARY- WAVE LEAKAGE 



The final result for the hexagon is then given in Eq. (k 



In this section we compute the leakage due to boundary 
waves. To estimate the loss we consider the total internal 
reflection {9i > 9c) of a plane wave at an infinitely ex- 
tended dielectric interface. This consideration is justified 
if Rc(fci?) ^ 1 and di not to close to the critical angle for 
total internal reflection 6c ■ Since 9e is a complex number 
in Eq. (|ll]), the wave in the lower-index region is evanes- 
cent, i.e. it decays exponentially with increasing distance 
from the boundary. Along the boundary the evanescent 
wave propagates with constant velocity. The correspond- 
ing total energy (or better intensity) flux at a given point 
at the boundary can be computed by means of 



Sxdy , 



(Bl) 



where Sx is the x-component of the Poynting vector |l2 



S = — Re(E X H*) = -ReMS/^P*) . (B2) 



Integration of Eq. (Bl) using Eqs. (B2) and ( |ll| ) yields 
c n sin 0, 1 



1 



1 



(B3) 



with a from Eq. (p^). 

Now we relate the total flux to the intensity inside the 
resonator. The intensity of the wave function tpi + ipr per 
area can be easily calculated to be l/Air. Assuming that 
the boundary waves fully leave the cavity at corners, the 
outgoing relative intensity per unit time is given by 



AJ, 



n sin 6j 



bw 



X 



^ k\J sin^ ( 



1 1 -I- 



(B4) 



where % is the number of corners (where the light is emit- 
ted), A is the area of the resonator covered by the family 
of long-lived rays. Formula (B4) is valid for any regu- 
lar polygon. For the hexagon, x — ^ ^J^d A is given by 
the area of the hexagon minus the region not accessible 
by the family of long-lived rays as depicted in Fig. |l^. 
Elementary geometry yields 



A = ^(i?2 



Rl) = 



(B5) 



APPENDIX C: PSEUDOINTEGRABLE 
LEAKAGE 



Here we compute the pseudointcgrable leakage. We 
put the angular shift A6' from Eq. (A5) into the ini- 
tial conditions for the ray dynamics. Elementary geom- 
etry, see Fig. Il^, shows that after one round trip, i.e. 6 




FIG. 17; Region (empty hexagon) not accessible by the family of 
long-lived rays. The region is bounded by a hexagonal "caustic" 
with side length _Rc = R/\/^. 



bounces, the angular deviation gives rise to the spatial 
deviation As = |2/A6'|. The time for each round trip is 
M = In/c. 

From Fig. |l^ it is clear that points on the boundary 
within the distance As from a corner leave the cavity 
after the next round trip. Hence, the relative outgoing 
intensity can be computed as 



1 As 
A/n = — T- 
P i? At 



Attc I3{n) 
3Re{k)R'^ rfi 



(CI) 
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